Evidence for TeV Gamma-Ray Emission from the Galactic Plane 
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Gamma-ray emission from a narrow band at the Galactic equator has previously been detected up 
to 30 GeV. We report evidence for a TeV gamma-ray signal from the Galactic plane by Milagro, a 
large field of view water Cherenkov detector for extensive air showers. An excess with a significance 
of 4.5 standard deviations has been observed from the region of Galactic longitude I G (40°, 100°) 
and latitude |6| < 5°. Under the assumption of a simple power law spectrum, with no cutoff, in the 
EGRET- Milagro energy range, the measured integral flux is cf> y (> 3.5TeV) = (6.4 ± 1.4±2.1) ■ 10~ n 
cm~ 2 s -1 sr -1 . This flux is consistent with an extrapolation of the EGRET spectrum between 1 and 
30 GeV in this Galactic region. 

PACS numbers: 95.30.Cq, 95.85.Pw, 96.40.Pq, 98.35.-a 



Gamma rays are the best direct probe of cosmic rays 
outside the solar neighborhood. The interstellar medium, 
with its relatively high density in the Galactic plane, acts 
as a passive target for gamma-ray production by ener- 
getic cosmic rays. Mechanisms include interactions with 
gas cloud nuclei that produce gamma rays via 7r° decay, 
as well as cosmic-ray electron bremsstrahlung and inverse 
Compton scattering with the interstellar radiation field. 
Emission from a diffuse source concentrated in the nar- 
row band along the Galactic equator was indeed detected 
by the space-borne detectors SAS 2, COS B and no- 
tably EGRET at energies up to 30 GeV. Above 1 
GeV, the EGRET data in the region of the Galactic cen- 
ter show a hard spectrum E~ a with a differential spec- 
tral index a w 2.3, and a flux enhancement of as much as 
60%, compared to models with tt° production as the sole 
mechanism Q and using the local cosmic ray spectrum. 
Models that predict a Galactic flux enhanced by up to an 
order of magnitude over the tt° mechanism at very high 
energies were proposed 0, IE IS 0- Upper limits have 
been set by several groups in the TeV range 0, 0, 
and above 180 TeV |ll|. 

The Milagro Gamma Ray Observatory fl2l fl3] | is a 

large-field-of-view telescope designed to detect gamma 
rays near 1 TeV using water Cherenkov techniques to 
observe air shower particles that survive to the ground 
level. It is located at a latitude of 36° and an altitude 
of 2630 m in the Jemez Mountains, New Mexico, USA. 



A 60m x 80m x 8m covered pond, filled with clear wa- 
ter, has a top layer of 450 photomultipliers (PMT), used 
to reconstruct the shower direction with angular resolu- 
tion of about 0.75° from the relative PMT timing. A 
bottom layer of 273 PMT's is used for discrimination 
between gamma-ray and the dominant hadron-induced 
air showers. A Compactness cut, described in |l3j . re- 
jects about 90% of the hadronic background and retains 
about 45% of the gamma-ray signal for typical gamma- 
ray source spectra. On the opposite side of this cut there 
are 9 times more cosmic rays and 1.2 times more gamma 
rays, so that any gamma signal is suppressed by a factor 
of 7.4, compared to the Compactness-cut signal. Addi- 
tional cuts in this analysis require a minimum of 50 top 
layer PMT hits, at least 20 of which participate in the 
angle fit, a zenith angle of < 50° and a declination of 
10° < DEC < 60°. Periods with abnormal event rate, 
zenith or azimuth distribution, corrupted records, or un- 
stable operation were excluded. The results presented 
are from three calendar years of data collected by the 
Milagro detector starting July 2000. 

Figure ^ shows the EGRET gamma-ray flux along 
the Galactic equator Q], peaking near the Galactic 
center, and the Milagro exposure. The Galactic cen- 
ter is not visible to Milagro. Superposed is the nom- 
inal expected relative significance, the product of the 
EGRET flux and the square root of the Milagro expo- 
sure. The region with highest expected significance, Rl 



2 



(/ E (40°, 100°)), and a second region of high Milagro 
exposure R2 (2 £ (140°, 200°)) both with \b\ < 5° were 
selected a priori, as two critical regions for separate sta- 
tistical tests of a TeV gamma-ray signal. The Milagro ex- 
posures in Rl and R2 are similar. The average EGRET 
flux in R2 is about a factor of two lower than in Rl and 
the extrapolation to TeV energies may well be different 
in Rl and R2. Rl includes the Cygnus arm and R2 is in 
the extreme outer galaxy. 
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FIG. 1: Top: 4 to 30 GeV EGRET flux, Bottom: Milagro 
Exposure (dashed), and expected relative significance (solid), 
along the Galactic equator. Regions Rl and R2 are indicated. 



Gamma-ray emission from a gamma-ray source should 
appear as an excess of the observed events (N s ) above the 
expected background count (Nb) in a candidate source 
direction bin. If the background is isotropic in the ce- 
lestial equatorial coordinate Right Ascension (RA) and 
the detector acceptance is not changing over some time 
window, the number of detected background events can 
be factorized into a time-independent acceptance shape 
G(x) in local angular coordinates x (e.g. Hour Angle 
and Declination) and a time- varying rate R(t). For a 
candidate source bin of interest (Rl, R2, or any bin of 
Figure 2) this gives N b = f ,. G(x) ■ R(t) dxdt. 

o / o u J source bin v ' v ' ■ ■ 

We take advantage of a "time-swapping" method [15j to 
perform a Monte Carlo calculation of this integral. The 
data is split into 8-hour segments in which the x and t of 
recorded events are considered random samples of G(x) 
and R(t). By pairing x and t randomly selected from 
these samples, new Galactic coordinates are generated 
for background events. Nb is incremented whenever such 
an event falls into the source bin. 

We use a self-consistent modification, described in [13, 
of the method outlined above. The assumption that G(x) 
is time independent is relaxed, allowing the incorpora- 
tion of a small diurnal modulation, from atmospheric 
changes, that is observed in the zenith distribution of 
events. The statistical error of the background is cal- 



culated with <JN b 2 = A ■ Nb + Nb/ B, where A is the 
mean on-source to off-source exposure ratio and B = 10 
is the number of time swaps per event |l7|. The events 
arriving from the Galactic plane are excluded from the 
"time-swapping" , in order to maintain the statistical in- 
dependence of N s and Nb ^3 ■ 

A very small anisotropy of the cosmic-ray background, 
at the level of a few parts in 10 4 , was observed in a sep- 
arate Milagro study |2J| and is applied as a correction 
in this analysis. The anisotropy is consistent in shape 
and magnitude with reports from several other experi- 
ments Milagro's observed anisotropy is de- 
scribed adequately by three longest wavelength harmon- 
ics in right ascension (RA), whose amplitude and phase 
vary linearly with declination (DEC). The anisotropy 
study provides a functional form for the background cor- 
rection in this analysis. The subtracted sky map in RA 
DEC coordinates, with Rl and R2 excluded, is fit to this 
form, then extrapolated into these regions. The resulting 
fractional correction of the background, applied in Table 
|H is S bg = -0.63±0.30-l(T 4 and S bg = ±0.04±0.30-lCT 4 
for Rl and R2 respectively; the errors are the statistical 
errors of the performed fit. 

The test for a source signal is done taking Rl and R2 as 
two separate single bins. The steps are given in Table [I] 
with the fractional excess defined by T = (N s — Nb)/Nb- 
The statistical errors on T and Sb g are combined in 
quadrature in the last line. An excess with a significance 
of 4.5 standard deviations is seen in Rl |25j, while no 
excess is detected in R2. 



Region 


Rl 


R2 


iV, 


238, 095, 657 ±15430 


254, 800,416 ±15962 


N„ 


238, 025, 840 ± 8003 


254, 826, 272 ± 8841 


N s -N b 


69, 817 ± 17, 382 


-25, 853 ± 18,247 


raw 


(2.93 ±0.73) ■ 10~ 4 


(-1.01 ±0.72) ■ 10~ 4 


Sbg 


(-0.63 ±0.30) ■ 10" 4 


(±0.04 ±0.30) • 10" 4 


T 


(3.56 ± 0.79) ■ 10~ 4 


(-1.05 ±0.78) ■ 10~ 4 



TABLE I: Calculation of the fractional excess T = (N 3 — 
Nt)/Nb for Rl and R2. N a is the observed event count; the 
error is on the mean of the underlying Poisson distribution. 

Figure[21(top) shows the two-dimensional color map of 
the significance in Galactic coordinates. An enhancement 
in a ridge along the Galactic equator in Rl is the most 
prominent feature. A Gaussian fit to the distribution 
of significances in Figure |2| (bottom) has the requisite 
normal distribution with zero mean and unit variance, 
while entries from Rl are shifted to the right. 

Profiles of the fractional excess T in latitude and longi- 
tude are shown in Figure|21 The peak of the enhancement 
is just north of the Galactic equator in Figures and 
Such a shift of the gamma signal towards positive lati- 
tudes was also seen in the EGRET data 0, interpreted 
as a large scale warp of the Galactic H 1 clouds |26| . 
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FIG. 2: Top (color online): Map with a color scale of the 
significance in Galactic coordinates in 5° x 5° bins. Bot- 
tom: Significance distribution outside RI (solid) and inside 
Rl (dashed). 
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FIG. 3: Profiles of the fractional excess in latitude for the 
Rl longitude band I £ (40°, 100°), and in longitude for the 
latitude band |6| < 5° of Rl and R2. A fit of the EGRET 
longitudinal source shape (solid line) superposed in Rl. 



To test for possible systematic problems in the analysis 
of the Rl excess, an identical analysis, including all mod- 
ifications and corrections, was performed on 14 months 
of data, with the Compactness cut reversed to impover- 
ish gamma rays relative to the cosmic-ray background. 
The result for Rl, T' = (0.68 ± 0.39) • 10~ 4 agrees well 



with ^suppressed = ^/7.4 = (0.48 ± 0.11) • 10" 4 , that is 
expected in the absence of systematic effects. 

We interpret the excess in Rl (Table HJ as dif- 
fuse gamma-ray emission, or emission from unresolved 
gamma-ray sources, or both. The ratio of integral 
gamma- and cosmic-ray fluxes is calculated with: 



(H+He) 



r](a 7 ) 



(1) 



Because the intrinsic energy resolution of Milagro is 
limited, individual event energies were not used in this 
analysis |27| . Our flux determination assumes that the 
gamma rays have a power law spectrum with no cutoff. 
The energy scale is determined from air-shower and de- 
tector simulations and confirmed by measurements of the 
Crab Nebula [13J . This gives a median energy of approxi- 
mately 3.5 TeV, with a 20% systematic error, for gamma 
rays from the Galactic equator. The coefficient 7]{a-y) is 
the energy and transit averaged ratio of gamma-ray to 
cosmic-ray effective area, with the Compactness cut, ob- 
tained with Monte-Carlo air-shower and detector simula- 
tions [i^ . We report results for the power law index con- 
necting the top EGRET point (10 - 30 GeV) and the Mi- 
lagro point at 3.5 TeV, a 7 = 2.61, for which 77 = 6.2±2.0. 
The error on r\ is the estimate of Monte-Carlo uncer- 
tainties, including the energy scale error. We note that 
the flux is only sensitive to the absolute energy scale of 
Milagro in proportion to the difference of power law in- 
dices E~( a ~t ~ acr '> . The cosmic-ray integral flux above 3.5 
TeV/nucleus is <j> H +He = 1.2 ■ VQ~ 6 cmT 2 s' 1 sr' 1 H<J. 

Integral flux results are shown in Table ITT1 and plotted 
in Figure H together with the EGRET data. 





SINGLE BIN 


MULTIBIN FIT 


K 


(5.7 ± 1.3 ± 1.8) ■ 10" s 


(5.7 ± 1.3 ± 1.8) • 10"° 




(6.8 ±1.5 ±2.2) ■ lO" 11 


(7.3 ±1.5 ±2.2) ■ 10" 11 


Q 7 


2.61 ±0.03 ±0.05 


2.60 ±0.03 ±0.05 



TABLE II: Ratio of the Milagro gamma-ray to cosmic ray 
flux, integral flux above 3.5 TeV (in cm~ 2 s -1 sr -1 ), and the 
connecting EGRET-Milagro power law index, for region Rl. 
The last column uses a fit of the longitude profile in Fig. [S] 
(bottom) to the EGRET source shape. The results are not 
sensitive to this choice. 



For comparison, we fit th e hig h end (1 GeV to 30 GeV) 
of the EGRET spectrum HEJ in Rl to a power law, ob- 
taining the index ci^^egret) = 2.51 ± 0.05, softer than 
the Galactic center value of a ss 2.3, quoted earlier. Ex- 
trapolations of the EGRET fits in Rl and R2, with their 
1 standard deviation error corridors, are superposed on 
Figure 0] 

As seen in Figure^! the emission from Rl at a 3.5 TeV 
median energy is consistent with the extrapolation from 
the high-end EGRET data between 1 and 30 GeV. The 
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FIG. 4: Integral flux results of Milagro and EGRET, with 99% 
c.l. upper limits from Whipple {I £ (38.5°, 41-5°), |6| < 2°) 
H|, HEGRA {I £ (38°, 43°), \b\ < 5°) [|, Tibet (I G (20°, 55°), 
|6| < 5°) E3. The display for region R2 is shifted by 10" 2 . 



absence of a detected excess in R2 requires a slight steep- 
ening of power law to a 1 (R2) > 2.66 (99% cl) compared 
to a 7 (i?l) = 2.61 ±0.03. The Milagro flux upper limit of 
r> 7 (> 3.5TeU) < 4 • lO^^m^s^sr 1 (99% cl) in R2 
reaches the EGRET extrapolation for that region. 

These findings do not require models predicting a 
strong enhancement of the diffuse flux compared to 
conventional mechanisms, such as an increased inverse 
Compton component |30|. a harder proton spectrum in 
the Galactic plane @,or contributions from unresolved 
supernova remnants |5j. Assuming a power law spec- 
trum for gamma rays with no cutoff below 10 TeV, we 
can exclude for Rl a hard spectrum with power law in- 
dex a 7 < 2.48 (99% cl). The results are consistent with 
a gamma ray power law index that asymptotically ap- 
proaches that of cosmic rays, as predicted if 7r° produc- 
tion becomes the main source of the gamma flux 

In summary, Milagro has observed a 4.5 standard de- 
viation excess in the mid-longitude galactic plane region 
Rl. The consistency of EGRET and Milagro data under 
the simple power law assumption reduces the motivation 
for speculation about more complicated gamma-ray spec- 
tra. With its one integral measurement, however, Mila- 
gro cannot rule out alternate models. For the many pos- 
sible multiparameter models, whose examination is be- 
yond the scope of this letter, the Milagro result provides 
one constraint in their parameter space. In a model with 
a simple continuous power law spectrum from EGRET 
to Milagro energies, this measurement is evidence for 
gamma-ray emission at TeV energies, with a flux that 
is consistent with an extrapolation from EGRET, and a 
median energy of about 3.5TeV for the detected gamma 



rays. 
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